FIELD OF THE INVENTION 

The present invention relates to measuring cardiac function intervals. 
BACKGROUND OF THE INVENTION 

It is known that alteration of the QT, QTc or RR interval on the electrocardiogram may be a 
marker for sudden death(^ "^^). Measurements of the QT interval are generally taken from a 
12-lead electrocardiogram where one to three heart beats are analyzed either individually or 
averaged(^^'^^). The 12-lead electrocardiogram provides only point-in-time data reflecting 
approximately 17 seconds of time that is required to inscribe a 12-lead ECG. The QT 
interval duration is dynamic, however, and can vary by upwards of 100 msec in a twenty four 
hour period(*^"^^). Thus the measurement of either a single or a few 12-lead ECGs sampled 
during 24 hours will miss the beat-to-beat dynamicity data that is inherent in the changes that 
occur. The dynamic data reflecting changes in ECG intervals is captured by longer 
recordings, generally 24 hours of continuous ECG data, referred to as either 24 hour 
ambulatory ECG (AECG) or Holter Monitoring('^'^^"^^). Heretofore, beat-to-beat ECG data, 
both short and long-term recordings, has been averaged due primarily to constraints in 
computing power. Unfortunately, averaging minimizes the understanding of the beat-to-beat 
variability inherent in QT interval data. Moreover, methods to analyze large data sets of 
cardiac intervals have been incomplete. For example the methods for beat-to-beat binning of 
QT and QTc intervals described by Callahan and Shell where limited to analysis of only 
outliers(^^), calculating the % of beats that exceed a certain threshold. The disclosures by 
Shell and Callahan do not teach a method to analyze central tendency, variance, kurtosis or 
other statistical properties of the histogram as appropriate for Gausian or non-Gaussian 
distributions. 

Increases in the QT and QTc interval measurements on a 12-lead Electrocardiogram (ECG) 
are associated with an increased risk of cardiac dysrhythmias and sudden cardiac death. See, 
for example, Algra(^^), Schwartz(^^'^^) and SawickiC"^^'"^^). The increased QTc interval length 
is associated with an increased risk of sudden death from all causes. The prolongation of the 
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QTc interval induced by pharmaceuticals has been associated with Torsade de Pointes and 
sudden death; the pharmaceutical induction of prolonged QTc intervals has formed the basis 
for removal of pharmaceuticals from the market. There is, however, no readily agreed upon 
method to measure the dynamic changes in the QTc interval, particularly for long term 
recordings of the ECG. 

While the resting 12-lead electrocardiogram may provide important spatial information 
regarding the status of ventricular repolarization, the use of a single 12-lead ECG measured 
randomly in time may disregard potentially important prognostic data regarding the 
dynamicity, temporal relationships, and circadian rhythms of the QT interval. 

It is known that the QT interval may undergo significant changes over both the short and 
long term due to circadian rhythms. See, for example, Yi, et alC*"^) who teach the association 
between circadian rhythm and sudden death associated with acute myocardial infarction. See 
also, for example, Callahan and Shell who describe a method to assess circadian changes in 
the QT interval. 

It is known that the QTc interval may undergo significant changes over both the shorter and 
longer term due to autonomic control. See, for example, Cappatto et al, Browne et al('^^), and 
Kautzner, et alC^''^^ demonstrated the relationship between sympathetic and vagal tone on 
the QT and QTc interval. 

Thus, a single 12-lead ECG taken at a given point in time may provide misleading and 
inaccurate cardiac risk data. Therefore, analysis of the QT interval for an entire 24-hour 
period, reflecting circadian and autonomic changes, may provide additional information 
regarding the risk of sudden death not available on the single, random 12-lead ECG. 

It is now possible to measure the QT interval on 24-hour Holter (AECG) recordingsC*^'^^'"^*'"^^" 

These measurements have generally been reported as averages over short time periods, 
typically between about 15 seconds and about five minutes, for example Molnar et al(^^'^^) or 
Yanaga, et al(^^). The use of averaged QT measurements may obscure significant short-term 
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variations in the QT intervals. Conversely, beat-to-beat measurements retain the natural 
variability data that may be important for calculating a patient's risk of dysrhythmia and 
sudden death. 

More recently beat-to-beat QT interval measurements have been used but methods to analyze 
the beat-to-beat changes have been incomplete. 

Although beat-to-beat variability of the QT interval has been described by Berger and others 

little is known regarding normal ranges in variability and measures of the QT interval 
over a 24-hour period using beat-to-beat measurements. 

Molnar and colleagues published a study that gives some indication of the dynamic range of 
the QT intervals using five minute averages and not beat-to-beat measurements(^^).. They 
reported a mean maximum QTc interval of 495 ms for normal subjects using 24-hour 
ambulatory monitoring. They also showed a mean intra-subject change of 95 ms. Molnar 
further reported six normal female subjects as having a maximum mean QTc interval 
measurement of more than 500 ms. These mean maximum measurements were taken over a 
five-minute period. 

The use of average QTc measurements obscures the dynamicity of individual beats. 
Measurements of central tendency, skewness and shape of histograms have not been used 
extensively to describe the relationship of QT and QTc measurements in histograms 
representing beat-to-beat QT, QTc or RR intervals. These measurements may be important to 
give an overall picture of the status of the subject. 

It is an objective of the present invention, in a preferred embodiment, to enable the 
assessment of the QT and QTc intervals and other cardiac function intervals on a beat-to-beat 
basis, providing a composite histogram of the individual beats with QT and QTc intervals. 

It is another objective of the present invention, in a preferred embodiment, to 
enable the measurement and assessment of the QT and QTc intervals and other 
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cardiac function intervals over an extended period of time, including not only 
periods of time greater than about one minute but also periods of time lasting 
at least 24 hours and even longer, in some cases. 

SUMMARY OF THE INVENTION 

In accordance with the present invention, in a preferred embodiment, this and other 
objectives are achieved by providing a method for analyzing beat-to-beat QT intervals from 
high-resolution Ambulatory Electrocardiographic monitoring (AECG) to detect the 
frequency distribution in a continuous AECG recording. Beat-to-beat QT and RR intervals 
may be measured to calculate beat-to-beat QTc. In a preferred embodiment, a composite of 
the entire frequency distribution of QT and QTc intervals taken from a set of observations 
with a common characteristic may be examined. Moreover, a composite of one characteristic 
may be statistically compared to a composite with other characteristics, including statistical 
methods that do not assume a normal distribution of the histogram. 

The present invention, in a preferred embodiment, provides a method to analyze beat-to-beat 
QT data, stratify the data according to a time-series bin-array, and create a composite of 
multiple histograms. This method and apparatus may be applicable to a wide variety of 
different subjects including, for example, normal subjects, subjects with the Inherited Long 
QT syndrome (ILQTS), and subjects exposed to drug titration . The statistical characteristics 
of a normal subject group can be compared to either a second group, or to individuals who 
have taken a drug, have potential congenital heart disease, have been exposed to an 
environmental toxin, or have a disease which could cause prolongation of the QTc interval 
such as diabetes mellitus. 

Further objects, advantages and other features of the present invention will be apparent to 
those skilled in the art upon reading the disclosure set forth herein. 
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DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

The following detailed description is of the best presently contemplated mode of carrying out 
the invention. This description is not to be taken in a limiting sense, but is made merely for 
the purpose of illustrating the general principles of the invention. The scope of the invention 
is best defined by the appended claims. 

In a preferred embodiment, standard 24-hour AECG recordings may be obtained using any 
commercially available Holter cassette tape recording device. An example of this type of 
device is a Reynolds Medical Tracker II (Reynolds Medical, Hertford UK) recorder. Also in 
the preferred embodiment, standard 24-hour AECG recordings can be obtained using 
commercially available digital Holter recorders that have a sufficient sample rate to allow 
detection and measurement of the cardiac intervals. For QT interval analysis, the sample rate 
can be between 128 and 2000 samples per seconds. The preferred embodiment is the use of a 
sample rate of at least of at least 1000 samples per second. These digital recorders must also 
be compatible with a Holter playback system that can produce beat-to beat interval 
measurements. An example of this type of device is the Reynolds Medical LifeCard CF 
recorder (Reynolds Medical, Hertford, UK). An example of the compatible Holter playback 
system is the Reynolds Medical Pathfinder 700 series. (Reynolds Medical, Hertford, UK). 
These recorders and playback systems are commercially available and need not be modified. 

Analog signals from the Holter cassette recordings may be digitized at 12-bit or higher 
resolution using a Holter playback system that has the ability to perform interval 
measurements. The Reynolds Medical Pathfinder 700 series Holter analyzer is an example 
of this type of equipment (Reynolds Medical, Hertford, UK). 

Using the digitized file of the electrocardiogram, a QT interval analysis may be accomplished 
in the following manner: The onset of a Q-wave (Qb) may be defined and a cursor may be 
placed at this point. The end of a T-wave (Te) may be defined and a second cursor may be 
placed at this point. The data from the digital file may then be replayed at 60-times normal 
time, while the cursors on the Qb and Te points may be monitored for stability. If either 
cursor wavers from the Qb or Te points, the cursors may be replaced and the affected portion 
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of the data may be reanalyzed. The QT interval may be defined as the time difference 
between the time points at Qb and Te. The QT intervals may be measured for the entire 
AECG recording on a beat-to-beat basis. Other analysis systems that display digital data 

I 

may be used. 

The peak of an R-wave may be detected and a third cursor may be placed (Rp). Accordingly, 
each QT interval may be matched with the preceding R-R interval. For a 24-hour recording, 
this may result in approximately 100,000 beats for which a QT interval and an R-R interval 
may be defined. The data may then be output to a high-speed computer for post-analysis 
processing. 

In the examples described herein use was made of AECG recordings from normal volunteers, 
subjects treated with placebo and subjects on-treatment in a drug treatment study, and 
recordings from subjects with inherited Long QT Syndrome (ILQT). These recordings help 
to demonstrate the potential effectiveness creating composite curves in accordance with the 
present invention. 

In the examples described herein QTc was calculated by removing a time-series of the QT 
and preceding R— R intervals to a high-speed computer with both a fast processor and 
adequate disk storage space. For each QT interval, a QTc may be calculated using a variety 
of correction factors for the QT interval including Bazett's correction formula, Fridericia 
correction formula and linear correction formula. 

The QT and QTc intervals may be individually placed in the bins according to their 
measurement as described in Shell and Callahan. In a preferred embodiment the composite 
curves are constructed by software programs that generate a time series of approximately 
100,000 data points long of RR/QT/ QTc triplets for each patient. Then the QTc data for 
each patient is binned in a histogram for that patient, finally, software is used to merge many 
patients' data into a composite data set (a "population") and to take means and standard 
deviations of this population (assuming normalcy of the data). Finally, more the data thus 
aggregated into two or more populations can then be compared, again using a combination of 
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software and procedures as described in Press et al(^^), against each other to check for 
statistical difference between these two or more populations. 

These aggregated population curves can then be used as a template for comparison against an 
single patient's binned histogram to determine what population (e.g. normal, inherited 
disorder, or drug-induced damaged) this particular patient belongs . The current embodiment 
assumes normal distributions, but this is not intrinsic to the method and more sophisticated 
distribution-distinguishing numerical analysis and statistics is declared here as well. The 
numeric procedures used are commonly described in Press, et al. 

EXAMPLE 1-NORMAL SUBJECT Figure 1. 

In this example, the 24 hour ambulatory ECG from a normal subject was analyzed. The 24- 
hour ECG was digitized. The QT and RR interval was determined for each beat using the 
Reynolds's analyzer. All extra beats were eliminated. All beats with prolonged QT intervals 
were inspected and artifact was eliminated. The QT and RR files were then used to construct 
a Histogram of QT and QTc intervals. The histogram was constructed with 10 msec intervals. 
The histogram of QTc intervals is depicted in figure l.The normal subject had a mean QT 
interval measurement of 358 msecs with a standard deviation of 37 msecs. the mean QTc 
measurement was 409 msecs with a standard deviation of 13 msecs. 

EXAMPLE 2 - Figure 2 

PROLONGED QTc Intervals in Inherited Long QT Syndrome. 

The Inherited Long QT Syndrome is a genetic defect of the heart's ion channels. The patients 
with Inherited Long QT Syndrome are known to have intermittent prolonged QTc intervals. 
Often, however, many of the heart beats of patients with inherited long QT syndrome are 
within the normal range and the identification of these patients cannot be made from a single 
conventional 12 lead ECG. Since these patients, often children, die suddenly, failure to detect 
the presence of the abnormal gene can lead to sudden death of the infant, child or young 
adult, an unnecessary death since treatment is available to prevent such sudden death. In this 
example., a child with a known gene defect underwent 24 hour ambulatory monitoring. The 
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ECG was digitized and the QT and RR intervals defined. The QT and QTc histograms are 
depicted in figure 2. The mean QTc was 450 msec with a standard deviation of 20 msec . 

EXAMPLE 3. -Figure 3 

PATIENTS WITH DRUG INDUCED LONG QT INTERVAL. 

Many drugs can prolong the QTc interval and the drug induced prolonged QTc interval is 
associated with an increased incidence of sudden death. Many drugs have been removed 
from the market because they prolong the QTc interval. Cisapride is a drug that can prolong 
the QTc interval. In figure 3, the QTc interval histogram is depicted in a patient taking 
cisapride. The mean QTc interval was 440 msec. 

EXAMPLE 4. - Figure 4 

COMPARISON OF THE QTc INTERVAL USING COMPOSITE CURVES IN NORMAL 
SUBJECTS TO THOSE WITH INHERITED LONG QT SYNDROME AND DRUG 
THERAPY 

In figure 4, a composite curve was generated from six normal subjects and a composite curve 
was generated from six subjects with known Inherited Long QT Syndrome. The mean QTc 
from the normal subjects was 409 msec +/-20 msec while the mean QTc from the patients 
with Inherited Long QT Syndrome was 475 msec +/- 35 msec and the drug therapy was 430 
msec +/- 0.40 msec. The kurtosis for the normals was 525 while the skewness was 767. The 
kurtosis for the patients with the gene defect was 1.24 while the skewness was 0.203. 

EXAMPLE 5 - Figure 5. 

COMPARISON OF THE QTc INTERVAL USING COMPOSITE CURVES IN A 
SUBJECT BEFORE AND AFTER DRUG INTERVENTION 

Since one of the important uses of this methodology is to compare the QTc interval before 
and after the use of a pharmaceutical that could prolong the QTc interval, we compared a 
group of patients before and after the administration of a pharmaceutical. The 19 patients 
had 24 hour ambulatory ECG monitoring before and after the administration of drug. The 
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composite curves before and after the administration of drug are depicted in figure 5. The 
mode before was 391 msec and after was 392 msec. Using a paired t-test the p- value was 
0.98. The total number of beats analyzed before treatment was 8.5 million beats and was 8.7 
million beats after treatment. The use of such composite curves generates large data sets that 
allow determination of difference/no difference in treatment sets with a high degree of 
statistical reliability. The conventional method to define differences would have analyzed 
between 50 and 3000 beats taken from 12 to 128 patients on resting 12-lead ECG. 

We then compared the means in these two groups using a standard T-test and analysis of 
variance. The p-value for the difference was less than 0.000008. The creation of the 
composite curves allows definitive differentiation of the three groups. 

EXAMPLE 6. - Figure 6 

COMPARISON OF A SINGLE INDIVIDUAL TO A COMPOSITE SET OF DATA 
Frequently, one is confi^onted with the problem of defining if a set of QTc data is derived 
from a set of normal data. In this example, a single individual with Inherited Long QT 
Syndrome was compared to a set of normal subjects (figure 6). In this example, the mean for 
the normal set was 408 msec while the mean for the patient with IQLT was 501 msec. Then 
the ILQT patient's histogram was compared to the normal set by use of either analysis of 
variance or Student t-test, the p-value was less than 0.00000001 indicating that the likelihood 
of the patient's histogram was sampled fi*om the same population set as the normal of less 
than one in a million. This degree of statistical reliability would form the basis of a 
diagnostic test for patients with suspected ILQT. In this case the composite curve was 
comprised of 533,354 beats compared to 94,996 beats for the individual patient histogram. If 
the patient had a 12-lead ECG, there would have been fewer than 20 beats available to 
compare the QTc interval to a mean normal that did not account for the beat-to-beat 
dynamicity of the QTc interval. This example shows how the composite curve invention can 
be used as a diagnostic test. 

In a preferred embodiment, the present invention represents a new method for quantifying the 
QT and QTc interval measurements over a period of time. The invention allows a 
quantitative comparison of two or more sets of QT or QTc intervals. For example, the 
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invention allows comparison of a group of patients before and after a drug. The method 
described allows application of a variety of statistical methods to define whether two or more 
sets of intervals are different from one another. 

In a preferred embodiment, the method and apparatus may make use of high-speed 
computer processors, and large capacity data-storage media. In a preferred embodiment a 
1 GHz Pentium IV processor with an 80-gigabyte hard drive may be used to analyze and 
store the large data files. Several custom-built software programs are used to generate a time 
series of approximately 100,000 data points of RR/QT/ QTc measurements for each patient. 
Then the QTc data for each patient is binned in a histogram for that particular patient, finally, 
a combination of software and procedures are used to merge many patients* data into a 
composite data set (a '^population") and to take means and standard deviations of this 
population (assuming normalcy of the data). Finally, the data thus aggregated into two or 
more populations can then be compared, again using a combination of custom software and 
procedures, against each other to check for statistical difference between these two or more 
populations. 

These aggregated population curves can then be used as a template for comparison against an 
single patient's binned histogram to determine what population (e.g. normal, inherited 
disorder, or drug-induced damaged) this particular patient belongs . The current embodiment 
assumes normal distributions, but this is not intrinsic to the method and more sophisticated 
distribution distinguishing numerical analysis and statistics is declared here as well. 

Composite QTc histogram measurements in accordance with the present invention allows for 
a quantitative assessment of the number of specified intervals, such as QT and QTc in a 24- 
hour AECG recording. 

The present invention, in a preferred embodiment, is directed to a method for the 
quantification of beat-to-beat QT and QTc interval measurements from ambulatory 
electrocardiographic recordings. 
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A QT binning technique in accordance with the present invention may be used to provide 
information about the effects of a pharmaceutical. For instance, in the example illustrated in 
FIGS 5, patients had two separate 24 hour AECG recordings. The first was a base line ECG 
recording. Then a pharmaceutical agent was was given to the patients in random order and 
the patients were monitored. Using a binning method and construction of composite curves 
in accordance with the present invention, an increase in the QT interval could be 
demonstrated better than by simply averaging or measuring a QT interval. 

Although the preferred embodiment of the present invention has been described herein with 
respect to measurement and analysis of the QT interval, it will be recognized that a method in 
accordance with the present invention may also be useful in the measurement and analysis of 
a wide variety of other ECG and related biologically significant intervals 

In a preferred embodiment, the method takes discreet measurements and discreet intervals 
and places them into a time series bin or an amplitude series bin. For example, all of the RR 
intervals in a sample could be selected and coded according to their length and then placed 
into bins. Each bin could be characterized by a frequency. The same analysis could be 
performed using any interval on the electrocardiogram. 

The presently disclosed embodiments are to be considered in all respects as illustrative and 
not restrictive, the scope of the invention being indicated by the appended claims, rather than 
the foregoing description, and all changes which come within the meaning and range of 
equivalency of the claims are therefore intended to be embraced therein. 
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